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In recent work we have shown that a black hole stacked on a brane escapes once it acquires a 
recoil velocity. This result was obtained in the probe-brane approximation, i.e., when the tension of 
the brane is negligibly small. Therefore, it is not clear whether the effect of the brane tension may 
prevent the black hole from escaping for small recoil velocities. The question is whether a critical 
escape velocity exists. Here, we analyze this problem by studying the interaction between a Dirac- 
Nambu-Goto brane and a black hole assuming adiabatic (quasi-static) evolution. By describing the 
brane in a fixed black hole spacetime, which restricts our conclusions to lowest order effects in the 
tension, we find that the critical escape velocity does not exist for co-dimension one branes, while 
it does for higher co-dimension branes. 



I. INTRODUCTION 

A striking prediction of brane world models with TeV 
Planck scale is that small black holes can be created in 
high energy collisions of particles at energies within the 
reach of forthcoming experiments 0, 0) S S IE Ifl Q- 
When two particles on a brane collide at a center of 
mass energy larger than the fundamental Planck scale 
with sufficiently small impact parameter, the system col- 
lapses and a black hole forms. Soon after formation, the 
black hole will start emitting Hawking radiation partly 
into lower dimensional fields localized on the brane, and 
partly into higher dimensional bulk modes [3] • The emis- 
sion of higher dimensional modes will cause the black hole 
to recoil into the extra dimensions if there is no symme- 
try that suppresses the recoil velocity, like Z2-symmetry 
in co-dimension one case. 

The previous simple considerations motivated us to 
study the interaction between a small black hole and 
a brane, paying particular attention to how the system 
evolves dynamically when some perturbation gives the 
black hole a velocity relative to the brane. This prob- 
lem was initially studied in Ref. Q , in the context of a 
scalar field model, and it was shown that the black hole 
is capable of escaping from the brane when it recoils due 
to emission of higher dimensional quanta, but no physi- 
cal mechanism for the escape was suggested. This result 
was obtained in the approximation that the tension of 
the brane is negligibly small. 

In Ref. Ic?! we have considered the same problem from 
the different perspective of studying the dynamics of 
Dirac-Nambu-Goto branes in black hole spacetimes (Rel- 
evant results concerning the static case in 4 dimensions 
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were studied in Ref. [ij). Our results confirmed the 
conclusion of Ref. Q and, in addition, suggested a mech- 
anism for the escape of the black hole based on the re- 
connection of the brane. Subsequently, in Ref. [IJ, we 
tested these results within a field theory model, where the 
brane was described by a domain wall in a scalar effective 
field theory. In this way we could take into account the 
recombination processes and fully illustrate the evolution 
and escape mechanism. All the previous results were ob- 
tained in the approximation that the tension of the brane 
has no effect. Although ignoring the tension is a reason- 
able assumption when the recoil velocity is large, it might 
not be so in the opposite case of small recoil velocity. 

In order to make the problem more definite, let us 
briefly repeat some relevant points discussed in Ref. [l3| ■ 
There, we investigated the dynamics of a brane in a fixed 
Schwarzschild spacetime background and found that the 
brane pinches and the black hole slides off the brane. 
Having completely ignored the self-gravity of the brane, 
the black hole continues to move without being pulled 
back. Once we take into account the effects of a finite 
brane tension, however, it is not clear if the pinching of 
the brane occurs before the black hole is pulled back. 
When the tension of the brane is large, this problem is 
more difficult to study since we need to take account of 
the deformation of the geometry caused by the gravity 
of the brane. Hence, we restrict our consideration to the 
effects which are lowest order in the brane tension next 
to the probe-brane approximation. 

Our strategy goes as follows. We start by consider- 
ing the quasi-static evolution of the brane, by which we 
mean that the brane evolves adiabatically through the 
sequence of static solutions shown in Fig. ^ In the next 
section, we will provide a justification for this. Once the 
sequence of solutions is fixed, it is possible to evaluate 
the energy for each configuration. This allows one to see 
how the energy changes along the path that describes the 
black hole escape, and to estimate the required initial ki- 
netic energy for the black hole to overcome the height of 
the energy barrier along the path. As one might notice. 



there are several non-trivial technical issues to address. 
The most important one concerns the evaluation of the 
energy of an infinitely extended object, which is diver- 
gent and thus requires regularization. This issue will be 
addressed in the next section, proposing a regularization 
scheme. The regularization scheme we adopt is to fix the 
circumferential radius of the brane at the outer boundary, 
which is let to infinity in the final steps of the calculation. 
Although it seems unnecessary to justify this regulariza- 
tion scheme since it is apparently the most natural way, 
it is not so easy to give a rigorous proof. In Appendix 
we give an argument to justify the proposed regulariza- 
tion scheme. In Sec. Illll we give an analytic comparison 
of the energies between the initial and the final configu- 
rations. In the initial configuration the brane is placed 
on the equatorial plane of the black hole, while in the 
final one the brane is disconnected and far away from 
the black hole. This analysis anticipates the shape of 
the potential along the escape path. In Sec. IIVI the en- 
ergy for the sequence of static configurations is evaluated 
numerically and, based on the shape of the potential, we 
discuss the existence of a critical escape velocity. We find 
that there is no barrier in the potential for co-dimension 
one branes. In other words, the force between a black 
hole and a co-dimension one brane is repulsive. Hence, 
the configuration with a higher dimensional black hole 
attached to such a co-dimension one brane is unstable. 
In contrast, when the co-dimension is equal to or greater 
than two, there is a potential barrier along the escape 
path of the black hole in the configuration space, and 
therefore there is a critical escape velocity. Namely, the 
configuration with a black hole on a brane is stable under 
a perturbation due to a small recoil velocity. 



II. ENERGY OF STATIC CONFIGURATIONS 

We consider a system which consists of a black hole 
and a brane. Assuming that the tension of the brane a is 
small, we neglect the gravitational perturbations caused 
by the presence of the brane. Namely, the spacetime 
is assumed to be given by d-dimensional Schwarzschild 
geometry [l^[l5| . 
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with /(r) — 1 ~ 1/r ^^, where we set the product of 
the d-dimensional Newton's constant Gd and the mass 
of the black hole m to unity. We treat the brane as a 
{p+ 1)— dimensional Dirac-Nambu-Goto membrane. Im- 
posing spherical symmetry on the brane, the action is 
given by 



S 



dtL 



a dt dr {r sin 9) 



p-i 



'l-—e^+r^f0'^ , 



where 9 is the azimuthal inclination angle, a = crEp-i 
and Ep_i is the area of a (p—1)— dimensional unit sphere. 



In the static case, the Lagrangian L is identified with the 
energy E as 



E = -L = a dr{r sin 6*)^"^ (l + r 
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It is convenient to describe both r and 9 as functions of 
a parameter A. Then the energy can be written as 



E= / dA£ , 



where 



£ = a(rsinf?f-^ ' ^"^'^ ^ ' -2 
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Using the re-parametrisation invariancc of the above ex- 
pression, we fixed the boundary values of the parameter 
A to be and 1, therefore the boundaries are located at 
{r(O),0(O)}and{r(l),0(l)}. 

Now we consider a sequence of static solutions, varying 
the location of the boundaries. By this variation, the 
energy changes as 



SE 



dZ 



■5r - 



d£ 



a{rsa\ 



■59 



d\) 
p-1 



(5) 



vW^^^vW 



dr 2fd6. 

dx^' + '^Tx' 



The previous equation relates the the difference in energy 
between two neighboring static configurations to the 'en- 
ergy flux' at the inner and outer boundaries. The bulk 
contribution to the difference in energy disappears with 
the aid of the equations of motion. 

Here we are interested in the sequence of solutions 
that describes the escape of black hole from the brane 
as shown in Fig. ^ For these configurations the inner 
boundary is then fixed on the horizon or on the axis at 
9 — 0, while the outer boundary is lifted, thus describing 
the recoil in the center of mass frame of the black hole. 

In estimating whether there is a potential barrier in 
the path of escape for the black hole, first we have to 
fix the sequence of configurations. Thus, given an ini- 
tial static configuration, we lift the outer boundary of an 
small amount and look for the next static configuration 
to which the brane will relax. Let us consider the equa- 
torial configuration as the initial one and lift the outer 
boundary, as shown in Fig. |21 We know that such an in- 
termediate configuration with lifted outer boundary does 
not exist as a stationary one. If we introduce a dissipation 
term to the dynamics, the brane will continue to relax to 
a configuration which realizes a smaller value of energy. 
By minimizing numerically the action with fixed inner 
and outer boundaries, one can see that the brane evolves 
towards the limiting piece-wise configuration presented 
in Fig. [21 which is the one with smallest action amongst 




FIG. 1: The picture shows the sequence of static configurations describing the escape path of the black hole with minimum 
energy. The plot is obtained by solving the equation of motion for the brane numerically in the static case. The results 
reported here refers to the cases with a d = 5, 6, 7, 8 dimensional bulk spacetime and a p = 3 brane. R and Z are the cylindrical 
coordinates defined in Q. 



all the configurations with the same boundary conditions. 
This configuration is composed of two pieces: one is the 
static solution which satisfies the regularity condition on 
the inner boundary and the other is the piece lying on 
the black hole horizon. To reach such a configuration 
it will take infinite time because of the small lapse near 
the horizon in Schwarzschild coordinates. In the quasi- 
static evolution, however, the lifting proceeds through a 
sequence of infinitesimal changes of location of the outer 
boundary. Therefore the piece-wise configuration will be 
achieved in a good approximation, except for the very 
vicinity of the horizon. The contribution to the energy 
from the piece near the horizon is negligibly small when 
the configuration is close to the limiting piece-wise one. 
(Notice that / = on the horizon, and dr/dX « near 
the horizon.) Now, the portion of the brane lying on the 
horizon does not contribute to the energy, therefore, as 
far as energy is concerned, the path for the escape will 
be described by the configurations plotted in Fig. ^ in 
which the location of the inner boundary is also varied. 



Let us now examine the change rate of the energy of 
the static configurations coming from the variation of 
the locations of the inner and outer boundaries. First 
we consider the inner boundary. If we take into account 
the piece lying on the horizon, we can say that the inner 
boundary is always fixed at the initial location. Therefore 
its contribution should trivially vanish. As mentioned 
above, we can consider the sequence of static solutions 
instead of the piece-wise limiting configurations. In this 
case the inner boundary is varied according to the change 
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FIG. 2: When the outer boundary is slightly lifted the initial 
configuration relaxed to the one with the least action which 
is the piece-wise limiting configuration drawn above. 



of the outer boundary. However, when the inner bound- 
ary is on the horizon, we have Svq — and /o = 0, where 
the subscript denotes the value calculated at A ~ 0. 
When the inner boundary is on the axis, 69 — and 
dr/dO = 0. Therefore according to the above expression 
the contribution to the difference in energy from the 
inner boundary vanishes: 5Eq — 0. 

More complicated is the evaluation of the contribution 
coming from the outer boundary, 5Ei . The total energy 
is divergent if we consider an infinitely extended object. 
As a regularization scheme, we cut the configurations at 
a finite radius. It is convenient to introduce cylindrical 
coordinates 

R^r sine, Z = rcos6l. (6) 

For later use, we record the expression of d6/dr in terms 



of dZ/dR: 
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dr 



l-tan^f 



' dZ 
KdR 



tan( 



(7) 



We consider the situation in which the outer boundary is 
varied with R fixed. Since R represents the circumferen- 
tial radius, fixing R provides a natural way to choose the 
boundary in a coordinate invariant manner. One could 
imagine an alternative cutoff radius by fixing the value of 
f^R with an arbitrary power v. Here we discuss only the 
simplest case with i/ = 0, but it is easy to check that any 
choice of the cutoff radius of this kind leads to the same 
conclusions. In this sense, the results obtained below are 
robust. In Appendix we give a little more rigorous jus- 
tification of this chosen prescription for a regularization 
scheme. 

Since the value of R is fixed at the outer boundary, 6ri 
and 56i are related to each other by the relation 



5Ri = smOiSri -f ri cos9i59i — 0, 



(8) 



where the subscript 1 denotes the value at A = 1. The 
variation of the energy, 5Ei , due to displacement of the 
outer boundary can be written as 



6Ei 



a{r sin( 
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'■fifrY 
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,de 

dr 



5Zi 



In the following, we assume that Z becomes almost 
constant at a large R. This is expected because the 
spacetime is asymptotically Minkowski. This assump- 
tion can be directly verified by numerical integration of 
the equation of motion, and it is visually clear already 
from Fig. ^ Under this assumption, we have dZ/dR ^ I 
and tan 9 — R/Z ^ f . Thus we find that the first term in 
parentheses in the denominator of the R.H.S. of Eq. {Tj) 
can be safely neglected. At a large distance we can make 
use of the Newtonian approximation. Namely, we can 
treat Sf as a small perturbation denoting / by I -|- 5f. 
Using these approximations, we have 



SEi 



R-1Z 



RP-' — 
dR 



SZi, 



(9) 



where q = d — {p + 1) represents the co-dimension of the 
brane. For any type of brane the first term vanishes in 
the limit i?i — > oo. Hence we have 



SEi 



i?l— >00 



*«rMS)/2. 



(10) 



If we assume that 8E\j8Z\ becomes a constant inde- 
pendent of R\ at a large distance, we obtain dZ/dR -^ 
a~^R^~P{5Ei/5Zi). Then we have the asymptotic be- 
haviour of Z as 



R'^-P SEi 
{p-2)dJz'i 



for J) > 3. Zac is the asymptotic value of Z. In the case 
with p = 2, Z continues to change logarithmically even 
at a large radius. Here we concentrate on the cases with 
P> 3. 

It is easy to see from numerical calculations and also 
it is possible to check analytically that the asymptotic 
form of the solution is indeed the one mentioned above. 
In fact, if we write down the equation of motion keeping 
the terms linear in Z up to 0{5f), we obtain 



d'^Z p-ldZ 



dR^ 



R dR 



= (?-!) 



Z 



R'i- 



(12) 



where we neglected the terms of 0{Sf) multiplied by the 
derivative of Z because Z is almost constant. The right 
hand side gives the leading order correction due to Sf. we 
can easily solve the above equation completely neglecting 
the right hand side as Z = Z^d + aR^~P , where Zoo and 
a are integration constants. This solution is completely 
consistent with the expression (|II|I . We can iteratively 
evaluate the effect due to the right hand side of Ea. (|I2|l . 
Substituting Z = Zoo on the right hand side, the cor- 
rection to Z is obtained as ?~\x ZooR'^~p~'^- Therefore 

q(d~3) "" 

this correction becomes sub-dominant compared with the 
term aR^~P for a sufficiently large R. Hence, the use of 
expression (|II|I is justified. 



III. ENERGY OF MINKOWSKI- TYPE BRANE 
AT A LARGE DISTANCE 

We now consider the energy of static Minkowski-type 
brane solutions. As we mentioned in the preceding sec- 
tion, we can use Newtonian approximation at large r. In 
other words we expand the solution Z{R) in terms of df 
as 

ZiR) = Z(o) + Z(i) + • • • . 

The lowest order solution is the one in Minkowski space- 
time and it is therefore given by Ziq-\ = constant. Zri\ is 
the correction of 0{6f). The expression for the energy 
E can be also expanded with respect to df formally as 



E — £'(0) + £'(1) 



where 



-S'(o) ~ ^ i dX{r sin ( 
Jo 




pl 2dl 



£(1) = - / dXirsm9)P-'- 

Then, we have 



-Jf. 



(II) E[Z] = £(0) [Z(o)] + (^J dX ^^(1)) + E^i) [Zio)] , 



where we have neglected terms of O (((5/)^). Since the 



second term vanishes because Z, 



(0) 



constant satisfies 



5E(^Q)/5Z = 0, we only need to evaluate i?(i)[Z(o)] as 
the lowest order correction of 0{Sf). Substituting Sf = 



1/r -^ and 



rcosi 



-(0) 



we have 



E(i)[Z^o)] = -^^('o7 I d0{cosey-'isinO^P-' 



_Ri^oo 



^aZ' 



, r[p/2]r[q/2] 

(0) 4r[{p + q)/2y 



(13) 



For any number of co-dimensions q > 1, the 0- integral 
converges in the limit i?i — > cxd. The behaviour of 
i?(i)[Z(o)] is peculiar for the cases with q = 1 and q — 2. 
In the case with g = 1, i<^(i)[Z(o)] continues to increase 
for an increasing value of ^(o)- In the case with q = 2, 
Eii-\[Z/Q\] converges to a constant value (j/2p. In other 
cases, £'(i)[2'(o)] vanishes in the limit Z(o) -^ oo. 

We can now compare the energy of this asymptotic 
configuration with the one on the equatorial plane. The 
difference i?(o)[^(o)] - i^cquatoriai is given by the area of 
the brane cut by the black hole horizon, d /p. Hence, 
neglecting the terms of 0{{Sf)^), we have 



V^r[p/2] 

-(^Zr^-rrr, ', ,,, tor g = 1 



'AT[{p+l)/2] 



2P 



for q = 2 , 
for q> 3 , 

for a large value of Z^o , where we have defined 
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AE = E[Zcx,\ — ^/equatorial 



(15) 



A way of interpreting the previous results is to write down 
the equation for the gravitational field: 



I — 'f^^u — ^ iif 



1 



rf-2 



9^iuT 



where T^^ is the stress tensor coming from the brane. By 
taking the 00 component we obtain an equation for the 
gravitational potential /iqo: 



D/i, 



00 



2{d-2) 



-T, 



00 



The coefficient in the previous expression describes what 
happens: for g = 1 the sign is negative and the potential 
is repulsive, for q > 2 the sign is positive and the potential 
is attractive. The case of g = 2 is marginal. 

These results are completely in harmony with those ob- 
tained in the opposite approximation in which the black 
hole is treated as a test particle on the background space- 
time determined by the gravity of the brane. In this 
approximation, co-dimension one branes cannot be em- 
bedded in an asymptotically flat spacetime. The sim- 
plest solution with the highest symmetry is the Vilenkin- 
Ipser-Sikivie model. We prepare two copies of Minkowski 



spacetime and consider the inside of the hyperboloid 
whose invariant distance from the origin is constant. 
Placing a brane on this hyperboloid means gluing to- 
gether two copies of Minkowski spacetime there. Since 
the bulk is simply given by Minkowski spacetime in this 
model, the trajectory of a test particle is just a straight 
line. On the other hand, the brane moves along a hyper- 
boloid, and hence the brane is always accelerated out- 
ward. To the contrary, if we look at a test particle moving 
in the bulk from the brane, it is accelerated in the direc- 
tion leaving from the brane. Namely, the force acting on 
the particle looks repulsive. 

When we consider a co-dimension two brane, we can re- 
call a cosmic string in four dimensional spacetime, which 
we know is described by Minkowski spacetime with a 
deficit angle. For co-dimension two branes in general di- 
mensions, the same is true. Therefore a particle moving 
far from the brane does not feel the gravity force due to 
the brane at all. This is consistent with the fact that the 
derivative of the energy d-B(i) [.Z(o)]/dZ(o) computed from 
Eq. H13|l vanishes for q = 2. For higher co-dimensions, 
gravity force caused by a brane is attractive and is pro- 
portional to £^~*, where £ is the distance between the 
brane and the particle. This dependence on £ also agrees 
with the estimate from the derivative of the energy H13|l. 



IV. POTENTIAL BARRIER FOR THE ESCAPE 
OF BLACK HOLE 

In the present section we will complement the analyt- 
ical results of the preceding section by numerical esti- 
mates. To start with, we can compute the profile of the 
potential along the path of escape of Fig.^by calculating 
the difference of the regularized energy for the sequence 
of static solutions numerically. Here, due to technical 
simplicity, instead of evaluating the energies of the solu- 
tions directly, we evaluate at finite radius the difference 
in energy for neighboring configurations by using formula 
®- Fig. O shows the plots of {Zoo, «} with 



1 SE 



a SZo 



(16) 



for various co-dimensions. The close-up views of Fig. ^ 
illustrate the transition point where the brane pinches 
and goes from a configuration intersecting the black hole 
to the detached one with Minkowski topology. Finally, 
in Fig. |S1 we plot the value of the energy AE vs Z^o ■ 
Here AE, defined by Eq. H15(l . is evaluated by the inte- 
gral AE = J dZoo{SE/SZoo)- The asymptotic values of 
AE in the limit Z^o ^^ oo are consistent with the ana- 
lytic estimates given in Eq. H14|) . Fig. |H1 illustrates the 
detailed structure of the profile for the energy in the re- 
gion near the pinching. Besides this tiny structure the 
curves are monotonically decreasing for co-dimension one 
branes and increasing for the other cases. The most im- 
portant point is that there is no energy barrier which is 
higher than the initial value for co-dimension one branes. 
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FIG. 3: The figure shows the plots of {Zoo, a} for a p = 3 brane in a d = 5, 6, 7, 8 dimensional bulk respectively. The red line 
corresponds to the brane intersecting the black hole, whereas the blue line corresponds to the brane outside the black hole. 
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FIG. 4: Transition between black hole topology to Minkowski topology. 



For the other cases with higher co-dimensions, the barrier 
height does not exceed the asymptotic value at Z^c — >■ oo. 
These features can be more easily shown by considering 
the energy of a momentarily static configuration with 
Z ^constant. The energy for this configuration is given 
by 



Ez=const. = a / dRRP-'Jl - Z^R^ + Z^) 



-{d-l)/2 



The plots of this simple expression are similar to but 
slightly bigger than the results shown in Fig. |S1 This 



simple expression gives an upper bound on the poten- 
tial obtained by correctly solving the equation of motion 
for static configurations. Therefore this simplified po- 
tential is sufficient to show that there is no significantly 
high bump in the profile of the potential. In Fig. |B| we 
can distinguish three branches describing the evolution of 
the brane: the (stable) black hole branch where the brane 
intersects the black hole horizon, the (stable) Minkowski 
branch where the brane lies outside the black hole horizon 
and finally an unstable branch joining the first two. This 
last unstable branch contains the configuration attached 
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FIG. 5: The figure shows the dependence of A_E on Zoo for d = 5, 6, 7, 8 and p ■ 




FIG. 6: Close up of the profile in the region near pinching. 
The red branch represents the brane configurations attached 
to the black hole; the blue branch represents the brane sepa- 
rated from the black hole. The plot refers to the case d = 6, 
p = 3. 




FIG. 7: Schematic representation of Fig. El 



to the pole at 6* = on the event horizon, which we call 
the pinching configuration. This configuration represents 
a kind of 'geometric threshold' between the configura- 
tions that intersect the black hole and those that do not. 
If we closely look at the vicinity of the pinching config- 
uration, there is a finer structure Q, which is beyond 
the resolution of the plot given in Fig. EJ Schematically, 
these three branches are drawn in Fig.[7| Accordingly, we 
can distinguish three energy scales: Epinching, Edetached 



p\d 


6 


7 


8 


9 


10 


11 


3 


0.1523 


0.2405 


0.2830 


0.3044 


0.3158 


0.3223 


4 


- 


0.1171 


0.1800 


0.2099 


0.2254 


0.2340 


5 


- 


- 


0.0962 


0.1477 


0.1706 


0.1822 


6 


- 


- 


- 


0.0820 


0.1251 


0.1436 


7 


- 


- 


- 


- 


0.0687 


0.1081 


8 


- 


- 


- 


- 


- 


0.0608 



and E, 



instability ■ 



Epinching labels the energy of the pinch- TABLE I: Values of E. 



ing configuration. The values of Epinching are tabulated 
in Table 1 for all dimensions and co-dimensions up to 
d= 11. 



Edetached represents the minimum energy of a configu- 
ration detached from the black hole. It is interesting 



pinching 



for d = Q, ..., 11 andp - 



to note that such 'detached' configuration cannot be 
reached through a sequence of static solutions starting 
from the initial equatorial configuration, but neverthe- 
less it marks the threshold for the separation between the 
brane of the black hole, which can in principle be reached 



8 



via some different (dynamical) deformation of the brane. 
Finally, Einstabiuty represents the value of the energy cor- 
responding to the upper cusp (see Fig. EJ), which marks 
an onset of instability. If the brane starts from the initial 
equatorial configuration and goes through the sequence 
of static solution it will reach such 'instability' point. The 
question is what happens when this instability point is 
passed over. If we introduce dissipation we can expect 
that the brane will evolve towards a configuration with 
lower energy lying on the Minkowski branch and there- 
fore a configuration separated from the black hole will 
be realized. As one can notice from Fig.|31 the difference 
between the values oi Epi„chtng, Edetached and Einstahiuty 
is very small. Some values of Edetached and Einstauuty 
are tabulated in Table 2 for p = 3 and d = 6, ..., 11. 



d 


6 


7 


8 


9 


10 


11 


^instability 


0.1574 


0.2491 


0.2926 


0.3145 


0.3264 


0.3333 


-t^detached 


0.1510 


0.2384 


0.2806 


0.3019 


0.3134 


0.3200 



TABLE II: Einstabiiity and Edetached for p = 3 and d 
6, ...,11 



As anticipated, we find that the co-dimension one 
brane is peculiar in the sense that it does not have a bar- 
rier for the escape of black hole. This means that there is 
no critical escape velocity in the present approximation 
where the brane tension is not very large. By contrast, 
all other types of branes have a barrier. Roughly speak- 
ing, the height of the energy barrier is 0(cr(G'dm)P/'^'^^^^) 
where we have explicitly written the d— dimensional New- 
ton constant Gd and the black hole mass m. Then, by 
equating the kinetic energy of the black hole rav^ jl with 
this energy barrier, the critical escape velocity v^ will be 

estimated as Vc ~ \J ^G^J, mf^^'J)/'''^^)^ Very inter- 

estingly, the critical escape velocity Ve is independent of 
the mass of the black hole for co-dimension two branes 
{q = 2). In the cases with higher co-dimensions Vc is 
smaller for a black hole with a larger mass, as is expected 
from the larger inertia of the black hole. 



V. DISCUSSION 

In this paper we investigated a system composed of a 
brane and a black hole when a perturbation gives the 
black hole a relative velocity with respect to the brane. 
The principal problem we are interested in is the phe- 
nomena of escape of the black hole. The main goal of our 
study is to clarify whether the effect of the brane tension 
may prevent the black hole from escaping for small recoil 
velocities and to understand whether a critical escape 
velocity exists. 

Here, we analyzed this problem by studying the inter- 
action between a Dirac-Nambu-Goto brane and a black 



hole assuming adiabatic (quasi-static) evolution. By 
'quasi-static' we mean that the system evolves very slowly 
going through configurations which are stationary. Tak- 
ing the brane lying on the equatorial plane of the black 
hole as initial configuration, we considered the brane in 
a fixed Schwarzschild spacetime background, neglecting 
the gravitational perturbation caused by the brane, ap- 
proximation which is adequate in the case of small ten- 
sion. Our strategy is to compute the energy for the se- 
quence of the static solutions describing the escape as 
shown in Fig. ^ directly and see how the energy changes 
along the easiest path for the black hole to escape. From 
this knowledge we can estimate in which cases there is 
a potential barrier that prevents the separation and the 
initial energy necessary to overcome such a barrier. 

We gave an analytic estimate for the height of the po- 
tential by comparing the energies of the initial and the 
final configurations, where 'final configuration' refers to 
the one with the black hole disconnected and being far 
away from the brane. We also computed the energy of 
each static configuration numerically, thus providing sup- 
port to the analytical estimates. From the calculated 
shape of the potential, we find that there is no barrier 
in the potential for co-dimension one branes, meaning 
that the configuration with a black hole attached to a co- 
dimension one brane is unstable. In contrast, when the 
co-dimension of the brane is equal to or greater than two, 
there is a potential barrier along the escape path of the 
black hole, and the critical escape velocity is evaluated 

to be 0(v'crGj;''"''"^"/^''~^^m(2-9)/(''-3)), where d and 
q are the spacetime dimension and brane co-dimension, 
respectively. 

In Ref. [13 , we studied the phenomena of escape treat- 
ing the brane as a domain wall in a scalar effective field 
theory. Since here we have focused on the case Dirac- 
Nambu-Goto branes, it is interesting to ask whether the 
results of this paper are also valid when the brane has 
a thickness. We expect that even in the case of thick 
branes described by field theoretical topological defects, 
the analytic comparison of the energies between the ini- 
tial and the final state will not change dramatically. For 
the final configuration the description by a Dirac-Nambu- 
Goto brane is always a good approximation. For the ini- 
tial configuration, the contribution to the energy from 
the region near the horizon may be corrected, especially 
when the thickness of the defect is not very small com- 
pared with the horizon radius. However, the height of 
the barrier in the intermediate stage is potentially quite 
different, and therefore a more rigorous study of the thick 
case is necessary to make a definite statement. 

As we mentioned in the introduction, one of the main 
motivations that led us to consider the problem of how a 
black hole interacts with a domain wall is related to the 
possibility of observing mini black holes in collider exper- 
iments or in cosmic rays. In this context Ref. [l3| points 
out that the escape of the black hole might provide a way 
of distinguishing between various models. Specifically, 
one possible difference might come out from the fact that. 



contrary to models with large extra dimensions 0, Q, 
warped models have additional Z2— symmetry Q, which 
makes the brane behave as if the tension were infinite, re- 
sulting in the impossibility for mini black holes to leave 
the brane. This is, however, not the case in models when 
the Z2— symmetry is relaxed, as Refs. 0, [13, uM have 
shown. The results presented here take a further step in 
this direction and suggest that also within models of the 
same class, without .Z2— symmetry, the co-dimensionality 
has a remarkable effect. If the co-dimension is 1, then 
even for small recoil velocity the black hole is expected 
to slide off the brane, contrary to the higher co-dimension 
case where there is an energy barrier that would prevent 
this from occurring. 

At this point it seems interesting to compare the recoil 
velocity due to Hawking radiation to the critical velocity. 
If we assume that the recoil is due to inhomogeneous 
Hawking emission, a rough estimate gives Vrec — 1/ v^, 
where N ~ m/Tji is the number of emitted particles 
and Th ~ r^ is the Hawking temperature. The horizon 

radius rn is given by (rn/M^"^)^^^ with Md being the 
higher dimensional Planck scale. Thus, we have 



'^recl ^c 



Ml' 



1/2 



(17) 



From this equation, we can compute the critical value for 
the mass at which the recoil velocity equals the critical 
velocity: 



M' 



rf-2 



fjp+i 



Thus we can say that, if the initial black hole mass is 
smaller than the critical mass, the recoil velocity is ex- 
pected to exceed the critical one before long, and there- 
fore the black hole will leave the brane during the early 
stage of its evaporation. On the other hand, if the mass 
of the black hole is larger than rricr, we expect that the 
black hole will escape from the brane in a subsequent 
stage after its mass becomes as small as rricr after losing 
a significant portion of its initial mass. 

Finally we want to mention the validity range of the 
probe-brane approximation. The probe-brane approxi- 
mation will require that the mass of the portion of the 
brane near the black hole horizon is much smaller than 
the black hole mass: 



ar^ <^ m . 



This gives a constraint 



TO > Md 



Mr 



This constraint is not so stringent. As far as the brane 
tension is small in the higher dimensional Planck units, 
our approximation will be kept to be valid throughout 
the above-mentioned process of black hole evaporation 
until its escape from the brane at m ~ rricr- 
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APPENDIX A: JUSTIFICATION OF THE 
REGULARIZATION PROCEDURE 

We present one argument for justification of using the 
cutoff at a fixed circumferential radius. One way to make 
the energy well defined is to fix the outer boundary con- 
figuration completely; not only the position of the brane 
but also the gravitational field. More precisely, for ex- 
ample, fixing the boundary induced metric will be ap- 
propriate. This is impossible in the present setup since 
the spacetime is exactly given by Schwarzschild space- 
time. If we fix the location of the brane at the boundary, 
a different static solution corresponds to a different posi- 
tion of the black hole. Then, the induced metric on the 
boundary changes. To fix the boundary metric, we need 
to include metric perturbations. 

For clarity, we can consider the Newtonian picture. We 
set the outer boundary aX r — Tq and the brane is on the 
equatorial plane on the boundary. Then, the perturba- 
tion of the Newtonian potential (j) necessary to compen- 
sate the change of the location of the central black hole 
will be evaluated as 



d-3 Z 



■5Z. 



(Al) 



Here Gdm is set to unity as before. The leading order 
perturbation of the gravitational potential will be given 
by a dipole field 54> = ^^iSr, which scales as -- l/^o"^ 

for large Tq . The total amount of correction due to this 
perturbation of the gravitational potential to the grav- 
itational binding energy of the brane will scale like 54> 
multiplied by the area of the brane ~ rg. Hence it scales 
as ~ l/j^o- We can therefore neglect this contribution in 
the limit ro — > 00. The correction to the energy due to 
gravitational perturbations is quadratic in 5(j). Hence, its 
contribution scales as ScjP' x r''^^ oc l/r^"^, and we can 
neglect it safely. When we take into account the gravity, 
we also need to consider gravitational potential induced 
by the brane. However, the background gravitational 
potential is chosen to be a vacuum solution. Then the 
leading order correction is quadratic. Therefore it scales 
like a^ . In the end, we claim that the effect of gravity can 
be completely neglected when the outer boundary radius 



10 



To is moved to infinity as far as we consider energy of 0{<t) 
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